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Abstract 

We study in this paper the finite-size effects of a non-periodic lattice on a lattice calcu- 
lation. To this end we use a finite lattice equipped with a central difference derivative 
with homogeneous boundary conditions to calculate the bosonic mass associated to the 
Schwinger model. We found that the homogeneous boundary conditions produce absence 
of fermion doubling and chiral invariance, but we also found that in the continuum limit 
this lattice model does not yield the correct value of the boson mass as other models 
do. We discuss the reasons for this and, as a result, the matrix which cause the fermion 
doubling problem is identified. 
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1 Introduction 



In addition to large computing resources, numerical lattice calculations require of a de- 
tailed analysis of finite lattice models to support the numerical calculations in order to 
have control over the finite-size effects. The results obtained by numerical calculations 
using finite lattice models of quantum-field problems may not be completely correct if 
these effects are not taken into account [TJ. 

In a previous paper [2] we have considered the Schwinger model on a finite and non- 
periodic lattice consisting of zeros of the Hermite polynomials that incorporates a well- 
behaved lattice derivative and a discrete Fourier transform. We found that if an infinite 
number of sites are taken into account, the mass of the boson mode calculated analytically 
by using this lattice becomes the correct one, i.e, ej yfjf, but it is about 40% less than this 
value if a finite lattice is used in the calculations. 

Good agreement of some finite-lattice calculations with continuum results has been re- 
ported elsewhere [31 HI [5J [6] . Besides the use of different techniques in these papers, it is 
also worth to note the use of periodic lattices or periodic boundary conditions, which are 
conditions usually assumed in lattice theory. 

Our aim in this letter is to study the finite-size effects of a non-periodic model on a lat- 
tice calculation. To this end we calculate the bosonic mass associated to the Schwinger 
model by using a naive lattice model consisting in a set of iV evenly spaced points and a 
standard central difference derivative with homogeneous boundary conditions. We choose 
homogeneous boundary conditions because the quantum fields are expected to vanish at 
±oo, whereas periodic boundary conditions are more natural in bounded manifolds. 

The use of homogeneous boundary conditions introduces important changes with re- 
spect to the lattice model with periodic boundary conditions. Chirality is maintained 
with no fermion doubling. The Nielsen-Ninomiya theorem [7] does not apply since the 
lattice discretization used here is not translationally invariant. However, the continuum 
value of the boson mass for the one-flavor Schwinger model can not be obtained. 



2 The standard derivative 

Let the lattice sites x n be defined in (-L, L) by x n = [— 1 + 2nj (N + l)]L, n = 1, . . . , N. 
The lattice spacing is a = 2L/(N + 1) and the boundary points are obviously, Xq = —L 
and £jv+i = L. The lattice length should be an increasing function of N. It is convenient 
to choose the dependence 




as we show below. Since a continuum problem is discretized only at the lattice sites, the 
homogeneous boundary conditions allow us to write the elements of the differentiation 
matrix D as 

Djk=ir( S i+i*-hw)> j,k = 1,2,..., N. (2) 
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Eventually, it will be necessary to have a diagonal form of this derivative. The usual 
Fourier transformation can not be used here, since we are not using periodic boundary 
conditions. Thus, the first thing to do is to obtain a diagonalizing matrix F for the 
derivative (jSJ). Because the Fourier transform diagonalizes the derivative, this matrix will 
be the discrete Fourier transform for this problem. It can be obtained by solving the 
eigenvalue problem DV k = ipkVk, k = 1, 2 . . . , N. Let v n {p k ) denote the nth component 
of the vector V k . Making the transformation u n (p k ) = i n v n (p k ), the eigenvalue problem 
for D converts into the finite recurrence equation 

U n+ i(p k ) + U n -l{pk) = PkUniPk)- 

This is the recurrence equation for the Tchebichef polynomials. Therefore u n (p k ) is pro- 
portional to one of these polynomials. The normalization constant can be found by using 
the Christoffel-Darboux formula [§]. Thus, the unitary matrix F that diagonalizes D and 
solves 

DV k = ip k V k (3) 

is given by 

W« V N+l \U N ^{p k )y 9-1,2,. ..,N, ( 4 ) 

where U n (£) is the nth Tchebichef polynomial of the second kind. The eigenvalues of — iD 
are the momenta p k = £, k /2a, where is the kth zero of Un(£) and therefore, they are 
not repeated. Thus, if we denote by P the diagonal matrix whose nonzero elements are 
Pi,P2, ■ ■ ■ ,Pn, we have that 

F\-iD)F = P. (5) 

It is important to notice that the lattice model used here is not translationally invariant 
and that the lattice momentum P does not lay on a torus. Thus, the Nielsen-Ninomiya 
theorem does not apply and this model does not exhibit fermion doubling. In fact, any 
pair of matrices (F, D) satisfying a one-to-one relation like (jSJ) gives a lattice model free 
from doublers and chirally invariant (if the lattice is not translationally invariant), which 
is also the case of the lattice studied in [S] (see also Sec 4. of [ID])- Below, in Sec. S]we 
show that the differentiation matrix which is related to the fermion doubling problem is 
a nearest-neighbor-like derivative A, which is actually a generator of forward-backward 
rotations and it is indeed diagonalized by the discrete Fourier transform 
Since F is unitary, the set of eigenvectors given in (j3J for k = 1,2, . . . , N, can be used to 
expand any other vector in its F-components. For N great enough and k fixed, the form 
of Q is _ 

where q, k = 1, 2 . . . , N. Note that q is the index associated to the position x q and k is 
the momentum index. A more useful form is 



? qk = I exp 



TT / k 1 

ix a — 

q a \N + 1 2 



— exp 



TT f k 1 

—ix a — 1 — 

9 a \N + 1 2 
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where x q = qa. In order to take into account the sign of Pk, the index k should be changed 
to k + (N + 1) /2, where k = -{N - l)/2, . . . , (N - 1) /2. This change yields 



1 



qk 



f2N 



exp 



2Xr 



7T k 

aN+1 



— exp 



a \ N + 1 



+ 1 



For a very large number of points x q located symmetrically about the origin, 



(7) 



(8) 



In our case, Ax g = a = 2L/(N + 1). Taking into account ([T]), we obtain that l/y/2N 
Ax q /27r and therefore, 



or in short, 







(9) 



where we have defined 

K = k/(N + 1), fc = — (iV - l)/2, . . . , (N — l)/2. 
This formula will be used below. 



3 The lattice model 

We take a lattice Schwinger model based on the Heisenberg equations and canonical 
quantization, such that the space variable is discretized and time is considered to remain 
continuous. Once the differentiation matrix D has been chosen and the diagonalizing 
matrix F has been found, the discretization of the Schwinger model can be done along 
the same lines as in [2]. As a consequence, many expressions involving D or F obtained 
in Sec. 2 of [2] remain the same in this case. We are considering a lattice consisting of 
a very large but finite number N of points. This means that the sums yielded by the 
discretization of the integrals are finite sums that can be written in matrix form. Thus, 
if the number of nodes N is an even integer, the differentiation matrix D is nonsingular 
and the lattice hamiltonian Hl can be written in the form 



H L = -p t D- 2 p + i&(cr z ®D)V, 
where p and \1/ are related by Gauss' law 



(10) 
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The rth value of the axial current j is given by 

and the equation for the charge density is 

{dlp) r - (D 2 p) r = -±[(Dj) r ,itD-'p]. (11) 

As in Ref. [2], it is necessary to compute the Schwinger commutator [jj, Pj>] to obtain the 
asymptotic form of the mass term. To obtain this commutator, we express it in terms 
of the F-transforms of jj and pf. We begin by expanding *$f q and p q = eW q ^ q in its 
F-components 

k 

p q = eJ2l a l a k'Fl q F g k> + b k b\,F qk Fl q }. 

k,k' 

Thus, the F-transform p k = J2 q ^IqPq becomes 

p k = eJ2 Wk^Fl, q F qk/ Fl q + b k// bi,F qk „Fi q Fi q }. (12) 

q,k',k" 

Now, we have to take into account the explicit form of the matrix F. By using ([6]) and 
a trigonometric identity is possible to reduce the triple products appearing in ([121 . For 
example we obtain that 

F k/ , q F qk 'F kq = —— [Fl„_ k , +kq + F A f„ +A ,,_ fcj(? — F k „ +k , +kA — F k „_ k ,_ kq j (13) 

and using and the fact that Ak ~ 1/N, 

/ ;F k ,i g F q k'F kg = —^-(S K " tK '~ K + 8 K "- K ' +K — 5 K "- K i- K — S K " iK ' +K (14) 
q 

Therefore, the F-transform p k of the charge density becomes 

pk = 2vr 2^ a «'-« afc ' + a_ K , +tt afc' - a_ R ,_ R a k > - a K , +K a k > 
k' 

- - a - K '+,-i a k' + (J_ K ,_ K _ x a y + a\ l+K _ x a k , (15) 
+ b K r- K b k , + b- K r +K b k , — b- K i^ K b\, — b K i +K b\, 

- 6«'- re -ifcJ., - fe-^+K-ifoJ./ + + b K > +K -ib\,), 

where we have integrated the delta functions of argument nn/a. The F-transform j k of 
the axial current can be obtained similarly. In order to calculate [j k , p k >], it is necessary to 
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write this commutator in terms of normal-ordered operators with respect to the vacuum 
state filled with particles of any negative energy and antiparticles of any positive energy, 
which is defined by 

|vac) = n4lKl°>- 

K<0 K>0 

The use of the anticommutation relations 

{4, a K /} = 5 KK >, {b[, b K >} = 5 KK >, (16) 

and a very lengthy calculation shows that the normal-ordered commutator [jk, p~k>] anni- 
hilates the vacuum state. In order to explain this, let us take only the first and third lines 
of (IT5|) (corresponding to the first line of delta functions of (TH1) and to the delta centered 
at the origin of (jUj)) and the analogous counterpart of jk- These terms yield a contribution 
to the normal-ordered form of \jk,Pk'\ acting on |vac) given by 

2 

6 (X + + j + + 

"o / X a K-m a -K'-m a -n-m a — n'—m a - K+m a - k' -m & K J rm Q'— k' — m 4" a _ K _ m a K'-m 
271 * — ' 

in 

— Qi _ n 1 — — n n — -I- n ^ — — 



(17) 



+ a \,+m a K'+m ~ b_ R , _ m b- K - m + & /e /_ m &-., t _ m ~~ ^L K / +m ^-K-m + & re ' +m ^-/c-m 
_i_ /it Z, _ ftt /, I J,t /, _ fct I, I J,t !, 

~ b K '- m b-K+m + b^_ K , +m b_ K+m — & re / +m &- re +m — b^_ K , _ m b K+m + b K ,_ m b K + m 
~ b_ K t +m b K + m +m^t+m)| vac )) 

where we have used the fact that (ae/27i) 2 Na/n = e 2 a/2ir. Since m takes negative and 
positive values, this expression vanish identically The same situation occurs when the 
full expression of pk and jk are taken into account in the complete calculation of the 
commutator yielding that 

L?fc,P fc /]|vac) = 0. 
Thus, the inverse F-transform is also zero 

\jk,Pk'} = 0, 

and therefore, the right-hand side of (TTTj) vanishes 

m 2 p r = ^[(Dj) r ,p t D- 2 p} = 0. (18) 

This means that the bosonic mass associated to the Schwinger model is zero in this 
lattice formulation. This result is a consequence of the index symmetry exhibited by (JHj), 
which comes from the fact that the diagonalizing matrix F for the derivative (j2j) satisfies 
relations like ( Fl3l) . which, finally, are due to the specific form of F and D. Note that this 
null result is independent of the number N of lattice sites and the same result is obtained 
in the continuum limit. We give an explanation of this in the next section. 



6 



4 Discussion 



Firstly, let us consider a. N x N matrix T> very similar to (j2J) 
by the elements Di 2 and D 2 \ 



which differs from D only 



V 



1 

2a 



/ \/2 

>/2 

1 

V : : 
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(19) 



This matrix can also be considered as a derivative for a problem with homogeneous 
boundary conditions if N is large. It is diagonalized by the unitary matrix which solves 
the eigenproblem Wk = iPk^k [H]- The elements of V are given by 



(V, 



k)q 



qk 



,-9-1 . 



2(1 -pg) T q ^(p k ) 
N(l + 5 k0 ) |TV-i(p fc )| ' 



1,2, 



(20) 



where p fc = i]k/2a, and % is the fcth zero of Nth Tchebichef polynomial of the first kind 
(77). The asymptotic form of ( 1201) is 



qk 



-i) 



N(l + 5 h 



cos • 



fed J 



(g - l)(fc- 1/2)7T 

N 



(21) 



g, fc = 1, 2, ■ • • ,N. The argument of the cosine function in (T2T1) takes the form qkn/(N + 
1) ~ qkn/N [the one of ([SJ)], if the zeros ^ of Un(0 are used instead of ^ to evaluate 

Let us consider now an infinite lattice, i.e., a lattice with an infinite number of nodes and 
boundaries at infinity. Let 5 be the central difference derivative [generalization of (J2J)]. 
Since there is not a first or last point in this lattice, we should write 5 without reference 
to the first (or last) indexes 



\ 



1 
-1 1 
0-10 



V 
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Any finite principal submatrix of 5 coincides with the principal submatrices formed by 
deleting the first row and column of D and T> for some N, which become identical. There- 
fore, the asymptotic form of F or the asymptotic form of J 7 , or both, can approximately 
diagonalize the matrix S, and a common choice is a linear combination of the sine and 
cosine functions [cf. Eqs. (EI) and ( f2~TT) ] to form the usual discrete Fourier transform, 
eliminating the symmetry which gives a zero mass term ([18]) for the Schwinger model in 
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this lattice formulation and recovering the correct value for this term [2]. Note that any 
finite matrix with a nonzero border can be confused with 5 in the asymptotic limit. This 
is the case of 



A = — 

2a 



/ 


1 


•• 


• 


Tl\ 


-1 





1 •• 


• 








-1 


•• 


• 











■• 


• 


1 







■• 


■ -1 





(22) 



where ± is discriminated according to the parity of N. This matrix is diagonalized by the 
discrete Fourier transform exp(—i2irjk/N)/^/N and their eigenvalues can be ordered as 



i 2itj 
-sm— , 
a N 



where the index j runs over integers or half-integers [TT] . This is the matrix which yields 
the fermion doubling since it bends the Brillouin zone onto itself. 

It is worth to notice that Eq. f|T3|) plays an important role in the calculation of the mass 
term and to remind that the matrix F is determined by the boundary conditions. For 
example, for the usual discrete Fourier transform, there is only one term in the right-hand 
side of (JT3J) . As another example, we can take the matrices (JF, D) instead (F,D) and 
calculate again the mass term. In this case, Eq. (fT51) becomes 



1 

2~N 



3~k"q3~ l^'^kq — <-« »t ( ■^"ik" — fe'4-fc n •^•"j-fc'—fc n ~l~ ri "I" 3~lrt 



rt 



rt 



k"-k'+k,q 



k"+k'-k,q 



k"+k'+k,q 



k"-k'-k,q 



and taking into account that the value m = should be excluded from the sums since 
the number of nodes N is an even number, Eq. (fTTj) changes to an expression that does 
not vanish identically: 

i^k'k $—k'k)/ 1 {fl^rrflrn ^— «— rrfl— n— m) {flm®"tn ^n—m^K—m) 
7T z — ' L 
m 

"i" iPm^m b— k— rrfl— k— m) (Pm®" 171 ^K—m^n—m) |vac). 

The a-terms of the sum applied to |vac) produce 

-k\vac), S ^2(a\ n a m - a\_ m a K - m )\vac) = fc|vac), 



Q'rrflm ^—K—rrfl— k— rn) | Vac) 



giving a zero contribution to the mass term. The same result is obtained for the 6-terms 
of the sum and the total contribution to the term mass is zero again. 
Let us remark finally that for a finite lattice only one diagonalizing matrix for the deriva- 
tive D (or T>) exists, and therefore, there is no way to get rid of the null result which 
becomes a finite-size effect. 
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